Flattening of Single-Particle Spectra in Strongly Correlated Electron Systems 
and the Violation of the Wiedemann-Pranz Law 
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The renormalization of the Wiedemann-Franz (WF) ratio in strongly correlated electron systems 
is analyzed within the Landau quasiparticle picture. We demonstrate that the WF law is violated: 
(i) at the quantum critical point, where the effective mass diverges, and (ii) beyond a point of 
fermion condensation, where the single-particle spectrum e{p) becomes flat. Results of the analysis 
are compared with available experimental data. 
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a. Introduction. In Landau theory,— a foundation 
of understanding of phenomena in condensed matter, 
Fermi liquid (FL) is treated as a system of interacting 
quasiparticles, having an ideal-Fermi-gas-hke momentum 
distribution, (hereafter we set the Boltzmann's constant 
kR = 1), 



n{p) = 1 + e 



(1) 



In the vicinity of the Fermi surface, the FL single-particle 
spectrum e(p), measured from the chemical potential /i, 
also has the ideal-Fermi-gas form 



e{p) = vf{p - Pf) = Pf{p - Pf)/M* 



(2) 



where the bare mass M is replaced by a T-independent 
effective mass M* . Therefore FL thermodynamic and 
transport properties differ from those of ideal Fermi gas 
merely by a numerical factor — a feature, inherent in 
atomic nuclei, liquid '^He and conventional metals. 

One of the most prominent justifications of the ap- 
plicability of the Landau-Migdal quasiparticle picture to 
metals is associated with the low-temperature behavior 
of the Lorenz number L{T), the ratio of the thermal con- 
ductivity k{T) to the product of temperature T and con- 
ductivity ct{T). The T — > limit Lq, derived by Som- 
merfeld well before the creation of FL theory, equals^i^ 



Ln = lim — - — ^ — = — TT 
T^Q a{T)T 3e2 



(3) 



Eq. ([3]), known as the Wiedemann-Franz (WF) 
law, holds in normal states of electron systems of 
metals^i^iii^i^ except for i) heavy-fermion metal a^^d^d^ 
CeNiSn and CeCoIns, ii) an electron-doped materiali^ 
Pr2-2:Cea;Cu04_y, and iii) an underdoped compoundi^ 
YbBa2Cu30j,. In CeNiSn, the experimental value of the 
reduced Lorenz number L{T)/Lo ~ 1.5 changes little at 
T < 1 K that rules out the relevance of phonons to the vi- 
olation of the WF law. In the electron-doped compound 
Pr2_2:Ce2:Cu04_j,, the departure of the ratio L{T)/Lq 
from 1 at T > 0.3 K is positive as well and even larger 
than in CeNiSn. 



Predictions of FL theory, including Eq. ([3]), fail in the 
vicinity of the so-called quantum critical point (QCP) 
where the effective mass M* diverges, since at the QCP, 
the FL spectrum ([J) becomes meaningless. In a stan- 
dard scenario of the QCP j^^d^ the divergence of the ef- 
fective mass is attributed to vanishing of the quasipar- 
ticle weight z in the single-particle state close to points 
of second-order phase transitions, implying that the FL 
quasiparticle picture of phenomena breaks down. In deal- 
ing with the WF law, a scenario of its violation, associ- 
ated with critical fluctuations, is recently advanced in 
Ref. [l3. However, the standard scenario of the QCP is 
flawedfi^ and therefore in this article, we employ a differ- 
ent topological scenario of the QCP, where the departure 
of the Lorenz number L{0) from the WF value is as- 
sociated with a rearrangement of single-particle degrees 
of freedom, a phenomenon described within the Landau 
quasiparticle picture. 

b. FL formulas for transport coefficients. Within 
Landau theory, relation between conductivities a and 
K, and the Seebeck thermoelectric coefficient S, has the 
form^> 

^^^^ 'SHT)^VM- (4) 



a(T)T 



Io{T) 



Here 



S{T) 



and 



4(T) 



T 



dv. 



(5) 



(6) 



where r is the collision time, dv is the volume element of 
momentum space, and n{p) is given by Eq. ([1]). 

Overwhelming contributions to the integrals Ik come 
from a narrow vicinity |e| ~ T of the Fermi surface. In 
Fermi liquids, obeying Landau theory, the Seebeck coef- 
ficient S{T) vanishes linearly with T at T ^ 0, similarly 
to the specific heat, given by 



C{T) 



.dn{p) 



(7) 
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Furthermore, in an ideal Fermi gas with impurities, 
the two quantities are connected with each other by 
relationi^i 



(8) 



where N^v is the Avogadro number. In interacting Fermi 
gases, the single-particle spectrum e{p) is not a parabolic 
function of p, and q(T) ^ —1, nevertheless, the propor- 
tionality between S{T) and C(T) holds, and therefore the 
contribution of the Seebeck coefficient to Eq. ([4]) remains 
minor. 

In conventional Fermi liquids, the group velocity can 
be factored out from the integrals ^ . The same is valid 
for the collision time r, depending at T ^ merely on 
impurity scattering. This yields'^'^ Ii{T = 0) = 0, and 
hiT 0)//o(r 0) = ttVS. Upon inserting these 
numbers into Eq. fl]) we do arrive at Eq. ([3]) that holds, 
even if several bands cross the Fermi surface simultane- 
ously. 

It is worth noting that top-quality samples, used in 
modern experiments, possess an extremely low residual 
resistivity pq < O.bfiflcm, and often at temperatures, at 
which measurements are carried out, electron-electron 
scattering comes into play, rendering the FL damping 
rate 7 ^ e-dependent^S 

7(e, T) = 7(0, 0) -f 7eeT2(l + {e/2nTf) , (9) 

that alters the Lorenz number L(T). To illustrate the 
change of its value we neglect the impurity scattering 
and set 7(0,0) = 0. In this case, the integrals Ik 



Ik 



1 



7eeT2 



e\^f de\ dn{e{p)) 



T 



dp 



dv 



de 



1 + (e/27rT)' 



(10) 

are calculated numerically, the results are inserted into 
Eq. (g]) to yield 



L{T)/Lo = 0.84 



(11) 



that agrees with results obtained in Refs. |21||22| by a dif- 
ferent method, giving rise to a slight decrease of L{T) 
at temperatures where the resistivity p(T) ceases to be 
constant. This suppression explains the departure of the 
experimental WF ratio from the FL value ([31), found in 
the heavy- fermion metal YbBa2Cu3 0j,. 

c. Violation of the WF law at the topological quan- 
tum critical point. In homogeneous Fermi liquids, the 
QCP was uncovered first in experimental studies of a 
two-dimensional (2D) electron gas i^'^'^'^i^^'^^ where the 
electron effective mass M*(p) diverges at a critical value 
Tc — 7.0 of a dimensionless parameter = \/2Me^ /pr- 
Results21 of microscopic calculations of 2D electron spec- 
tra e{p, T = 0, rs) at = 5, 6 and 7 are shown in Fig. [1] 
We see that e{p, T = 0, r^) becomes flatter and flatter, as 
the critical value = 7.0 approaches, at which for the 
first time the group velocity v = de{p)/dp vanishes at the 




FIG. 1: Single-particle spectrum e(p) of a homogeneous 2D 
electron gas in units of e% = p%/2M, evaluated^^ at T = 
for different values of the dimensionless parameter Vs = 
V2Me^/pF. 



Fermi surface. As seen from Fig. [U the QCP electron 
spectrum, denoted further ec{p,T — 0), has an inflection 
point £cip, T = 0) (X (p — pf)^- 

At finite T, the QCP group velocity dec{p,T)/dp ac- 
quires a finite value,— and the QCP single-particle spec- 
trum becomes^- 



^c{p,T) =PF 



P-PF 

M*{T) 



V2 



+ ^ip-PFr 



(12) 



The T-dependence of the effective mass M*{T) is evalu- 
ated on the base of the Landau equatioiJ^ 



dejp) 
dp 



p If 



dn{pi) pjdpi 

dpi TT^ 



(13) 



connecting e{p) with the quasiparticle momentum dis- 
tribution n(j)), given by Eq. ([1]), in terms of the first 
harmonic /i of the interaction function /. If /(g) is an 
analytical function of q, leading T-dependent contribu- 
tions to e(p, T) come from first terms of Taylor expansion 
of f{p,pi) vs (jp — pf) and [pi — pp)- As a result, one 
obtains^S 



M 



M*{T) 



'TM 



V Pf 



2/3 



(14) 



To evaluate the integrals Ik , we introduce a new integra- 
tion variable e(p) instead of p, then express p — pp in 
terms of energy e from Eq. to find 




where y = p — Pf and yo — [iMpp /v2Y^^ ■ With these 
results, the integrals Ik are calculated numerically that 
yields 



iQCp(0)/Lo = 1.81 . 



(16) 



Thus close to the QCP, the ratio Lqcp (0) /Lq is enhanced 
as compared with the FL value ([3]). This result is in 
agreement with the experimental value of the violation 
of the WF law, observed in Refs. [l^El 
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d. Numerical calculations of the Lorenz number L(T) 
through the topological QCP. In this paragraph we show 
results of numerical calculations of Eq. (fT3|) with an in- 
teraction function /, having the analytical form 



A 



1 



(17) 



where the parameter A is fixed, while the parameter a, 
depending on r^, is chosen to provide the best agreement 
with the microscopic spectrum, drawn in Fig.[TJ Results 
of calculations of the momentum distribution n{p, T) and 
spectrum e{p, T) at different values of a are given in the 
second and third columns of Fig. [5] respectively. On the 
FL side, (the upper panels (a) and (b)), of the QCP, lo- 
cated ata = 0.479, the spectrum e(p, T) has the standard 
FL form ^ that holds until the effective mass M* at- 
tains values ~ lO^M. The QCP spectrum tc{p,T), with 
the group velocity vf{T) vanishing at T = is shown in 
the third panel (c). It differs from the FL one in two as- 
pects. First, at T ^ it has the inflection point. Second, 
its form drastically changes with T elevation. 

Having at hand the spectrum e(p), the transport in- 
tegrals © are straightforwardly evaluated. Results 
of calculations of the functions L{T) and K{T) = 
3/2(r)/(7r2/o(r)) are shown in the first column of Fig. H 
The upper panel (a) of Fig. [5] illustrates the situation on 
the FL side of the QCP (a = 0.520, M*/M = 6). We 
observe no deviations from the FL predictions. As seen 
from the panel (b), the WF law also holds in the im- 
mediate vicinity of the QCP (a = 0A82,M*/M = 80). 
However, the departure from FL theory, associated with 
the temperature dependence of the spectrum e(p, T), be- 
comes well pronounced in both the ratios K — 3I2/{tt'^Io) 
and S = Ii/Iq already at extremely low T ~ lO^'^e^. 
The QCP results are shown in the middle panel (c). 
We see that the WF law is, indeed, violated: the ra- 
tio L{Q)/Lq turns out to be in excess of 1, in agreement 
with the above result pB]) . 

Results, shown in two lower panels (d) and (e), where 
a = 0.475 and a = 0.472 correspondingly, demonstrate 
that as the system goes away from the QCP, the value of 
the group velocity vf{T = Q) turns out to be finite again, 
and the WF law ([3]) is recovered. However, already at 
T ~ lQ~'^e^p, even lower, than on the FL side, the Lorenz 
number L[T) becomes T— dependent. 

e. Topological phase transitions in strongly correlated 
Fermi systems. As seen from the panels (d) and (e), 
beyond the QCP, equation 



e(p) = l^ 



(18) 



has three roots pi < P2 < Ps, the curve e(j>,T — 0) 
crosses the Fermi level three times, and occupation num- 
bers n{p,T — 0) become: n{p) = 1 at p < pi, while at 
Pi < p < p2, n{p) = 0; at p2 < p < Ps, once again 
n{p) = 1, and at p > p3, nip) — 0. Thus at the coupling 
constant g > 5qcp, the Fermi surface becomes multi- 
connected. This is a typical topological phase transi- 
tion, at which no one symmetry, inherent in the ground 
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FIG. 2: The reduced Lorenz number L(T)/Lo together with 
the ratio K{T) = ?,h[T) / {-k"^ h{T)) (left panels) as functions 
of the reduced temperature T/e%, the occupation numbers 
n(p) (middle panels) and the single-particle spectrum (right 
panels) calculated at different line-type-coded temperatures 
with the interaction function (|17|) for five different sets of 
input parameters, which correspond to five different cases: 
Fermi liquid with the effective mass M*/M ~ 6, i.e. quite 
far from the QCP (panels a); Fermi liquid with M*/M ~ 80, 
which is close to the QCP (panels b); system at the QCP 
(panels c), the state just beyond the QCP, with P2 — pi — 
0.14pF (panels d) and the state well beyond the QCP, with 
P2 — Pi — 0.31 pf (panels e). 



state, is violated i^'^'^'^ As the interaction strength 
increases, the number of the points, where the spectrum 
e{p, T ~ 0) crosses the Fermi level, rapidly grows, how- 
ever, the number of roots of Eq. (|18p remains countable. 

In another type of the topological transitions, the so- 
called fermion condensation )^"^i^^i'^^i^^i'^^i^^ the roots of 
Eq. (fT5|) form an uncountable set. Indeed, the ground 
state energy i5 is a functionalii^ of the quasiparticle mo- 
mentum distribution n(p), confined within the interval 
< n < 1. In the strong coupling limit, the minimum of 
this functional is found from the variational condition"^- 



SE 
5n{p) 



= peC 



(19) 



the chemical potential /x, being determined by the re- 
quirement '^n{p) = p. Since the l.h.s. of Eq. p9p is 
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nothing but the quasiparticle energy e{p), we see that in 
the case at issue, solutions of Eq. (HH]), called the fermion 
condensate (FC), exist in a whole domain C, giving rise to 
swelling of the Fermi surface. Remarkably, the presence 
of the FC results in breaking of the particle-hole symme- 
try, another salient feature of the phenomenon of fermion 
condensation, that exhibits itself in a marked violation of 
the WF law. 

True, if the interaction strength is small, the solutions 
of Eq. ((T9| do not meet the restriction n{p) < 1. How- 
ever, beginning with a critical constant, i.e. at g > 5fCj 
these solutions, smooth hmctions n^(p), meet this re- 
striction wherever. Thus on the Lifshitz phase diagram, 
the FL phase occupies a region g < gqcp, a domain 
5QCP < 9 < 9fc is occupied by the phase, having the 
multi-connected Fermi surface, while a domain g > gpCj 
by the FC phase. 



follows we address a model with the interaction function 
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FIG. 3: Same as in Fig. [2] calculated "with the interaction 
function (|20[) for Ppp = 3 and four values of the parameter 
g in units of M^^: 3.0 (Fermi liquid, panels a), 3.55 (QCP, 
panels b), 3.9 (the state "with the FC fraction Pfc/p = 0-1) 
and 4.8 (the state with ppc/P ~ 0-5, panels d). 



f{q) = ge^p{-f3q)/q 



(20) 



The structure of the FC, emerging at g > gqcp, is 
found with the help of formulas, given in Ref. 30. It is 
shown in Figs. [3] and U) As seen from Figs. [3] and [H the 
model derivative dn{p)/dp has peaks at both the bound- 
ary points of the FC region. Notably vicinities of these 
points contribute overwhelmingly to the transport inte- 
grals Ik, since the single-particle spectrum e{p,T — 0) 
identically vanishes inside the FC domain. Outside this 
domain J2. 



de{p.T = Q) 
dp 



(21) 



This power behavior results in a marked violation of the 
WF law. Indeed, upon inserting Eq. (PTjl into Eq. one 
finds that major contributions to every of the integrals 
Ik, proportional to T^/^, come from the exterior of the 
FC domain, while its interior ensures minor ones, pro- 
portional to T . This conclusion is confirmed by results 
of numerical calculations, shown in Fig. [5] 
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/. Breaking of the particle-hole symmetry and the WF 
law in systems with a fermion condensate. It is worth 
noting that if the interaction function f{q) possesses a 
singularity at g = 0, then the Lifshitz phase diagram, 
constructed above, alters, since the phase with the multi- 
connected Fermi surface disappears, and at 5 > .gqcp 
the system contains the normal quasiparticles and a FC 
fraction that grows linearly with the difference g — gqcp- 
This situation is convenient for the demonstration of the 
impact of breaking of the particle-hole symmetry in sys- 
tems with a FC on the violation of the WF law. In what 



FIG. 4: Same as in Fig. [3] calculated for j3pF = 30 and four 
values of the parameter g in units of M~^: 25.0 (Fermi liquid, 
panels a), 33.4 (QCP, panels b), 69.0 (the state with p-pc/ p = 
0.1, panels c), and 224.0 (the state with ppcl P = 0.5, panels 
d). 

The transport integrals Ik are calculated numerically. 
In Fig. [5l we show results of these calculations, per- 
formed with the value of the dimensionless parameters 
[3pF — 3 and j3pF = 30. As seen from this figure, at 
the QCP and beyond it the integrals /o(T = 0) and 
/i(T = 0) has the same order, implying that in contrast 
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FIG. 5: The transport integrals Jq (solid lines), Ji (dashed 
lines) and I2 (short-dashed lines) in log-log scale as functions 
of reduced temperature T/e%, calculated "with the interaction 
function (|20|) for jSpp ~ S (left column) and /3pF = 30 (right 
column) and four values of the parameter g corresponding to 
FL (upper panels), the QCP (second line of panels), and to 
the states with 10% (third line) and 50% (lo"wer panels) of 
quasiparticles in the FC. 

to FL theory, the T = value of the Seebeck coeffi- 
cient S{0) — /i(0)/e/o(0) differs from 0, an exhibition of 
breaking of the particle-hole symmetry, occurring in any 
system with a FC. On the other hand, the ratio L(0)/Lo 
turns out to be even larger than at the QCP. Further- 
more, calculations demonstrate that with increasing /?, 
the value of this ratio increases as well, i.e. the longer 
the radius of the interaction function (|20|) in the coor- 
dinate space, the larger is the departure from the WF 
law. 

g. Anisotropy of the violation of the WF law close 
to the QCP in heavy-fermion metals. The anisotropic 
Fermi surface of the majority of heavy-fermion metals has 
a sector, where at T ^ 0, the quasiparticle group velocity 
vp keeps a finite value, implying that the WF law holds. 
However, recently in experimental studies of CeCoIn5 in 
external magnetic fields, close to the critical value He, 



suppressing superconductivity of this metal, the WF law 
was found to be violated,— The violation is anisotropic 
that cannot be attributed to the collapse of collective 
degrees of freedom. On the other hand, close to the 
topological QCP, the conductivity tensors (x <ViVk> 
and Kik oc <e{p)viVk> become anisotropic, and this 
anisotropy is well pronounced in sufficiently large mag- 
netic fields. Indeed, the magnetic field does not affect 
the z— components of the velocity v, parallel to its direc- 
tion. As a result, the particular QCP T-dependence of 
the integrals Ik holds, triggering the violation of the WF 
relation L^z = fzz/ (Tk^z) — 7r^/3e^. On the other hand, 
the electron motion is completely rearranged in the di- 
rection, perpendicular to H that leads to a considerable 
increase of the respective components of the group veloc- 
ity and the suppression of the departure of the respective 
components of the ratio L from their WF value. 

h. Conclusion. In conclusion, we have demon- 
strated that flattening of single-particle spectra e(p) 
of strongly correlated electron systems considerably 
changes their transport properties, especially beyond the 
point of fermion condensation due to breaking of the 
particle-hole symmetry. Results of our analysis demon- 
strate that search for the violation of the WF law in new 
materials should be confined to electron systems with- 
out disorder, possessing, nevertheless, a sufficiently large 
resistivity, the fact, justifying a minor role of light car- 
riers. One more distinguished feature of systems, where 
marked departures from the WF law can exist, is related 
to the enhancement of the Seebeck coefficient S{T 0) 
and the ratio q{T ^ 0) = eS{T Q)/C{T 0). Inter- 
estingly, these features are inherent in the heavy-fermion 
metals where the violation of the WF law was observed. 
Indeed, in CeNiSn, the q value exceeds 10^.^^ Its value is 
markedly enhanced in compounds Pr2-2;Ce2;Cu04_j^ as 
welliS 
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